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I. INTRODUCTION 

Gravitational radiation reaction is a problem of great 
interest in the detection of gravitational waves. For sec- 
ond and third generations of gravitational wave detec- 
tors, a leading candidate source is the radiation-reaction 
induced inspiral and merger of two compact objects like 
black holes or neutron stars. Moreover, the effects of 
spins are important for the emission of gravitational 
waves from such systems. Thus in order to develop highly 
accurate theoretical templates for gravitational wave de- 
tectors, one must study the gravitational radiation reac- 
tion from compact binary systems with spin effects. 

In the present paper, the leading-order post-Newtonian 
(PN) spin-orbit and spin(l)-spin(2) radiation-reaction, 
i.e., dissipative, Hamiltonians are calculated. This is the 
continuation of previous work in [1], where the formal- 
ism was prepared, and also extends the calculation of the 
3.5PN point-mass Hamiltonian in [2] to that of spinning 
objects. The contributions of the spin-dependent Hamil- 
tonians derived in the present paper to the equations of 
motion are 2.5PN orders weaker than the corresponding 
leading-order conservative ones. Recently, the contribu- 
tions to the motion of spinning objects have just come 
within reach of experimental verifications [3, 4]. A fur- 
ther increase in precision of experimental tests of general 
relativity will become available by creation and subse- 
quent improvement of gravitational wave astronomy [5, 6] 
in the future. For compact binary systems detectable by 
gravitational wave detectors, the Hamiltonians derived 
in the present paper become relevant in the late inspi- 
ral phase if one or more of the binary's constituents is 
rapidly rotating. And rapidly rotating black holes have 
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been proved to be astrophysically realistic [7]. There- 
fore, the derivation of the Hamiltonians with spin effects 
is necessary for the detection of gravitational waves. 

In this paper, we utilize the canonical formalism of 
Arnowitt, Deser, and Misner (ADM), which has not 
only shown to be valuable for calculating the conser- 
vative dynamics within the post-Newtonian and post- 
Minkowskian approximations (see, e.g., [8-10]) but also 
for the dissipative part of the dynamics [2] (with mis- 
prints corrected in [11]). Notice that the ADM formal- 
ism was extended from point-masses to objects with spins 
only recently [12] (see also [1, 13, 14]). This extension 
is valid to linear order in the single spins of the ob- 
jects, which not only includes spin-orbit but also spin(l)- 
spin(2) interactions. The remarkable structure of the ex- 
tended ADM formalism of the inclusion of the matter into 
the canonical field momentum [see Eq. (2.6)] is passing 
an excellent test in the present paper. For Hamiltonians 
of higher orders in spins see [13, 15-18]. 

Energy and angular momentum flux relevant for the 
PN order in question has been well known (see [19], for 
the next-to-leading-order calculation see [20, 21]). Based 
on these results, secular equations of motion for the or- 
bital elements corresponding to the leading-order spin- 
orbit and spin(l)-spin(2) radiation-reaction equations of 
motion were obtained in [22-24]. The general equations 
of motion at this order were calculated in [25-27] within 
the harmonic gauge. The Hamiltonians calculated in the 
present paper provide a compact expression which con- 
tains these general equations of motion (but within a dif- 
ferent gauge). And most importantly, the results in the 
present paper are valid for arbitrary many object sys- 
tems. The derived Hamiltonians are then applied to the 
calculation of the energy loss of a binary system, which 
is then compared with the well-known energy flux as a 
check. 

The conservative leading-order (PN) spin interactions 
for self-gravitating objects were derived some time ago 
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[28-31], see also [32, 33]. For the leading-order spin(l)- 
spin(l) radiation-reaction level calculations see, e.g., 
[23, 33, 34]. However, only recently the conservative 
next-to-leading-order spin effects could be treated, start- 
ing with the spin-orbit equations of motion in harmonic 
gauge [35] (with some extensions and misprints cor- 
rected in [36]). A corresponding conservative Hamil- 
tonian in the ADM gauge was obtained in [37]. The 
complete next-to-leading-order spin(l)-spin(2) conserva- 
tive Hamiltonian was hrst given in [38]. Other deriva- 
tions of the conservative next-to-leading-order spin-orbit 
and spin(l)-spin(2) dynamics can be found in [39-43] and 
a generalization to arbitrary many objects succeeded in 
[44]. Notice that the results given in the present paper 
are already valid for arbitrary many objects. Also the 
conservative next-to-leading-order spin(l)-spin(l) inter- 
action of black hole and/or neutron star binaries was de- 
rived recently [16-18, 45, 46]. The latter requires a mod- 
eling of the spin-induced quadrupole deformation, sec 
[33, 47]. Very recently, the conservative spin-dependent 
part of the post-Newtonian Hamiltonian was extended 
even to next-to-next-to-leading order for both the spin- 
orbit [48] and the spin(l)-spin(2) [49] cases. A poten- 
tial for the spin(l)-spin(2) case was simultaneously calcu- 
lated within an effective field theory approach [50]. No- 
tice that the conservative next-to-next-to-leading-order 
spin(l)-spin(2) Hamiltonian and the spin-orbit radiation- 
reaction Hamiltonian derived in the present paper are 
both of the order 4PN for maximally rotating objects. 
However, not all spin-dependent Hamiltonians up to 4PN 
for maximally rotating objects are known yet. We will 
in most cases use the phrase "formal n-th PN order" to 
represent our counting of PN orders in the present paper. 
This gives PN orders different from the maximally rotat- 
ing case, which we also occasionally refer to in the present 
paper (for a more detailed discussion see, e.g., Appendix 
A of [1]). But one should be aware that the spins are 
in fact further (independent) expansion variables. Spin 
effects were also considered within the post-Minkowskian 
approximation [51, 52]. 

The paper is organized as follows. First the ADM for- 
malism is reviewed in Sec. II. Then formal expressions 
for the radiation-reaction Hamiltonians in question are 
derived in Sec. III. Integrals appearing in these formal 
expressions are performed in Sec. IV. In Sec. V, the de- 
rived Hamiltonians are applied to the calculation of the 
energy loss, which is then compared with the energy flux. 
Finally, conclusions are given in Sec. VI. 

Our units are such that c = 1, but for the Newtonian 
gravitational constant G no convention will be used. This 
allows an easy transition to the different conventions for 
G used in [2] and [1] . For the signature of spacetime, we 
choose +2. Latin indices from the beginning of the al- 
phabet, such as a, b, label the individual objects. Greek 
indices run over 0, 1, 2, 3. Latin indices from the middle of 
the alphabet run over 1,2,3. Round brackets around an 
index denote a local basis, while round brackets around 
a number denote the formal order in c _1 , as in [1, 14]. 



A 3-vector x l is also denoted by x. Square brackets de- 
note index antisymmetrization and round brackets index 
symmetrization, i.e., a^b^ = \(a^b u + a v b^). 



II. THE ADM FORMALISM 

In this section, we provide a short overview of the 
ADM canonical formalism after gauge fixing [53] , see also 
[54, 55]. The Hamiltonian is given by the ADM energy 
expressed in terms of certain canonical variables, which 
also requires a (at least approximate) solution of the field 
constraints. 

The constraints of the gravitational field read 
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with the definitions 



^matter 

(2.1) 
(2.2) 



7r« = -^7(7* V - f 3 l M )K kl , (2.3) 

^matter = ^ffT^rfrf , ( 2 .4) 
^matter = _^ T .^ . ( 2 .5) 

They arise as certain projections of the Einstein field 
equations with respect to a timelike unit 4-vector 
with components n M = (-N, 0, 0, 0) or = (1, -N l )/N. 
Here, 7^ is the induced three-dimensional metric of the 
hypersurfaces orthogonal to n^, 7 its determinant, R 
the three-dimensional Ricci scalar, = —(57^,0 — 

Nu.j\)/N the extrinsic curvature, N the lapse function, 
N l the shift vector, the stress-energy tensor den- 

sity of the matter system, and semicolon denotes the 
three-dimensional covariant derivative. Partial coordi- 
nate derivatives di are also indicated by commas. 

For nonspinning objects, X ^ G ^ is the canonical mo- 
mentum conjugate to 7^ before gauge fixing. For spin- 
ning objects, the canonical field momentum has to be 
adapted, see [1, 12]. We write 

''"can = n 3 + ''matter ' (2-6) 

where ^matter contains spin-corrections. Throughout this 
paper we use the ADM transverse-traceless (TT) gauge, 
which is defined by: 



(2.7) 

(2.8) 



Here, 6{j is the Kronecker delta. And one has the decom- 
positions: 



7y ; = ( 1 + i ) S ij + h i? > 



can 'can 



I " rat 



(2.9) 
(2.10) 
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Notice that i, j, and k, etc., run over 1, 2, 3, and upper or 
lower an index is from now on done with the flat metric, 
thus changes nothing in the equations. We will ignore 
the difference of upper and lower indexes later and any 
two identical indexes can contract with no need to be 



one upper and one lower 
traceless, e.g., hJ, T 



hJ^ T and 7r*^ T are transverse- 



h TT - 



vector potentials V^ an and 



"ly "can 

= 0, and TT l J a 
*can b y : 



is related to the 



III. RADIATION-REACTION HAMILTONIANS 
UP TO FORMAL 3.5PN LEVEL 

In this section, we generalize the derivation of the 
radiation- reaction Hamiltonians up to the formal 3.5PN 
level performed in [2] so it becomes applicable to the 
spinning case. 
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It holds that: 
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(2.12) 



(2.13) 

(2.14) 
(2.15) 



A. Interaction Hamiltonian and wave equation 

We split the ADM Hamiltonian 7?adm into matter, 
field, and interaction parts, i.e., 



ii rrmattc] 
-"ADM — tl 

where the matter part H" 



H 



field 



H 1 ' 



(3.1) 



r is independent of the 
(truly dynamical) canonical field variables hj^ and 
7r^ T , the field part 7J ficld is independent of the canonical 
matter variables and reads explicitly: 
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field 
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(^) 2 + i^Tf 



(3.2) 



with the inverse Laplacian A 1 and 



5™ = §[(*« - A-'dAKSjk - A~ l djdk) 



{6 ik - A~ 1 didk)(Sjt - A-^di) (2.16) 
($ki - A~ 1 d k di)(S i . 



A^dtdj) 



The canonical field variables after gauge fixing are hj^ 
and ngl T . 

In order to obtain the ADM Hamiltonian, the four field 
constraints must be solved for the four variables </> and 
7f* an in terms of hJ 3 T , 7r*{^ T and canonical matter vari- 
ables, which enter through the stress-energy tensor via 
the source terms % mattcr and ^f or ^ ne K near or _ 

der source terms in spin, see [1, 12, 13] and also [14]). The 
canonical matter variables are the canonical position 
momentum P a j, and spin-tensor S a u\rj) of the a-th ob- 
ject. An analytic solution for <fr and 7r* an , however, can 
in general only be given in some approximation scheme. 
The ADM Hamiltonian is then given by: 
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ADM 



where S a ^) = ^t iik S, 



and 



'-ijk 



(2.17) 
is the completely 



antisymmetric Lcvi-Civita symbol. Hadm is the ADM 
energy expressed in terms of the canonical variables men- 
tioned above. The Poisson brackets read 

{^ T (x), **£V)} = lftrG^«*(x - x') , (2.18) 



{S a (i), Sa(j)} — ^ijkSa(k) , 

all others are zero. 



(2.19) 
(2.20) 



and the interaction part H mt depends on both canonical 
matter and field variables. The interaction Hamiltonian 
up to and including the formal 3.5PN level reads ([1], see 
also [2]) 



H 



int 
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d 6 x 
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%w ) ^ T 



+ 2 (^(3)^(2),j " 1"(5) ro atter) 7r c a r J 



(3.3) 



where 



B(4)y = 167rG 



b ( / d 3 xH^ s f tcv ] 



J 1 



<f>W,i<f>(2),j , (3-4) 



and V? 3 \ is a field quantity which will be discussed in 

Sec. IV. B/Q-jij is given by a similar expression [see 
(5.14) in [1] 



For comparison with [2], notice that 

J ki V(3)vw; - (0(2)^(3)) [Tffj) is another 

field quantity which will be discussed later]. Further, 



2B(4)jj is used in this paper. 



in [2] the quantity Au^j 

The equations of motion for the canonical field vari- 
ables follow from the ADM Hamiltonian by virtue of the 
Poisson brackets (2.18) as: 

1 LTT _ rTTij $H ABM -s 

16irG can ~ kl 5hJ? • [ ' 

Here the dot over a variable denotes the partial time 
derivative dt = J^. For quantities not depending on the 
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hypersurface coordinate x, this is to be understood as 
the ordinary time derivative. In terms of the interaction 
Hamiltonian H lnt , the field equations read 



1 nfrTT _ rTTij 

WttG ^ ~ dkl 



5H b 



d 5H' mt 

at 5tt^ t 
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r ijTT 



-—h TT 
16ttG 13 



A-/ 



can 

(5iJ int 

17 "can 



(3.7) 
(3.8) 



with □ = A — &1- To arrive at these expressions, the 
explicit form of H Rcld is used as in Eq. (3.2). Notice 
that it is easier to implement the boundary condition of 
no incoming gravitational radiation for a wave equation 
like Eq. (3.7) than for a system of first-order differential 
equations like Eqs. (3.5) and (3.6). Inserting the 3.5PN- 
accurate interaction Hamiltonian [Eq. (3.3)], one gets 



UhJ7 = 5 T J 13 



IB, 



(4)fc/ 



2B 



o —-/~iu matter lTT 
(6)kl - OTTUtip) <i kl 



MTT _ ±fcTT 
can rj ij 



+ 0(2) V™) 
f 



k-l 



(3) ( P(2),i - 7r (5)matter 



with the definition: 



B 



(6)ij 



(5.) matter 



(3.9) 
(3.10) 

(3.11) 



One can get alternative expressions for -B(4)ij and 
B(R\ij in terms of 7y = <J*jTij by comparing the wave 
equation for /i^ T with the Einstein equations (see [1]), 
e.g., 

= -87rG7(4)jj - -0( 2 ),i0(2) j ■ (3-12) 
This should agree with Eq. (3.4) after the TT-projection. 

B. Near-zone expansion 

At the considered order, aspect like tail effects play no 
role (see e.g. [56]). We may therefore solve the wave equa- 
tion for hj? by an order-by-order evaluation of the re- 
tarded solution. Further, the field solution is only needed 
in the near-zone. 

In order to discuss the near-zone expansion, we write 
the wave equation for hJ T schematically as: 



kl 



(3.13) 



The near-zone expansion of the retarded solution to this 
equation corresponds to a series in c _1 entering through 
the retarded time t ret = t — c _1 |x — x'|, reading 



LoSki — LiSki 
L%Ski — LsSu + ■ ■ ■ 



(3.14) 



where the TT-projector was pulled in front of the re- 
tarded solution and the integral operator L n is defined 
by 

(L n /)(x,t) = -.J-j |dV|x-x'rV(x',i). (3.15) 

Notice that L 2n = A -1- ™ for n G N, in particular Lq = 
A- 1 . 

Using the PN-expanded source of the wave equation 
from Eq. (3.9) one may arrange the near-zone expansion 
by PN orders as: 



_ , TT , iTT , h TT , h TT 

n ij - n (4)ij + n (5)ij + n (6)ij + n (7)ij 



(3.16) 



It is important that only a finite number of terms from 
the near-zone expansion [Eq. (3.14)] contribute to a spe- 
cific PN order due to the increasing number of time 
derivatives therein. Therefore, one obtains 



TT 



(3.17) 
(3.18) 

hj 7 ^ = tl Uj + tl 2ij + II.;, , + 'tiuj + Qij , (3.19) 



h(5)ij — X(4)ij i 
I 

where 



n 



lij 



26™ j L lB 



{4)kl , 
(6)fci j 



U 2 ij = 8n 



r 1 A TT ^ T fh TT O/ matter 
UO kl ^1 y n (4)kl H {2) 



n 3 y = 26™* Li {y^ m 



n. 



matter 
(2) 



(3.20) 
(3.21) 
(3.22) 

(3.23) 
(3.24) 
(3.25) 



Notice that the application of L\ to a total divergence like 
(</)( 2 )/i T ; T m ) im leads to a vanishing result. It will become 
apparent in the next section that hj^j is not needed in 
the present paper (but it contributes to the conservative 
3PN Hamiltonian). The definitions Puj = tluj, P 2 ij = 
^2ij, Puj = n 3 ij, and Rij = ti^ were used in [2]. 

An application of the operator L\ obviously leads to 
a field depending on time only (i.e., not depending on 
x). This allows an easy calculation of the (regularized) 
TT-projections in Eqs. (3.20) - (3.23) by means of the 
formula 



(3.26) 



valid for an arbitrary x-independent function Aki (t) (see 
[2]). Here STF denotes the symmetric trace-free part, 
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4™ = ^(A ij + A ji )-^S ij A kk . 



(3.27) 



Further, hJJ- is a function of time only, /iXT. ■ i = 0. As 
a consequence of these simplifications, we finally have 

*(*)« = ( 3 - 28 ) 
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Qij 



j3„, oSTF 
5^ / d ^(ejy- 



4G 



j3„ lTT nj matter 



(3.29) 
(3.30) 

1 /• / \ STF 

-^J d 3 x{V^ 2lj ) , (3.31) 
^5™' / dV|x-x'| 2 B (4)fc; (x',t), (3.32) 
^uC^W , (3.33) 



where the PN-expanded Hamilton constraint in the form 
^0(2) — — l6TrGH"^ ttcl was used to arrive at the last 
equation. 



C. Radiation-reaction Hamiltonians 

The dissipation through emission of gravitational radi- 
ation enters the PN-expansion via hj^j and hj^- , which 
are antisymmetric under time reversal. The parts of 
the Hamiltonian linear in hj^j or hj^j thus give the 
radiation-reaction Hamiltonians at the considered order. 
Notice that H Rcld does not contribute to the matter equa- 
tions of motion, so we only need to consider H lnt . The 
radiation-reaction Hamiltonians are thus given by: 



rrint 
-"2.5PN 
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H. 



3.5PN 



16ttG 
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16nG 



d 3 x 



(3.34) 



B, 



(6)ij 



/f^^aymattor i,TT \ iTT 



d 



16ttG dt 



j3„lTT tj 

a x «(5)ij7r( 5 ) matter ■ 



(3.35) 



where we used hJZ j>k = 0, with Eqs. (3.11) and (3.10). 
Equation (3.10) reads explicitly: 



(3.36) 



The last term in Eq. (3.35) corresponds to a canonical 
transformation and could be dropped, but we keep it for 
now. 

One has to be aware of a subtlety here. The matter 
variables entering the Hamiltonian via the solution for 
hj? play a special role as they may not be treated as 
dynamical (i.e., phase space) variables. Otherwise, the 
matter equations of motion resulting from the Hamil- 
tonian would in general be wrong (at the conservative 
level one can use a Routhian to avoid this problem, 
see [8]). Instead these nondynamical matter variables 
entering through hj^ are treated as functions depend- 
ing explicitly on time only. This introduces an explicit 
time-dependence into the radiation-reaction Hamiltoni- 
ans, which is a very natural description of a dissipative 
system via canonical methods. 



In order to distinguish the nondynamical matter vari- 
ables from the dynamical ones, we attach a prime to 
their object label as in, e.g., Pi' or P a /, and also talk 
of primed and unprimcd variables for short. Further, we 
introduce an explicit time derivative dj x , which only acts 
on the primed variables (The partial and ordinary time 
derivatives act on both primed and unprimed variables 
here). A superscript a — >• a' is attached to a field to de- 
note that its solution should be expressed in terms of the 
primed variables. This denotes an exchange of all object 
labels by labels with a prime, not just of label a. Thus 
hj^ - and hj^ - in Eqs. (3.34) and (3.35) should better 

be denoted by hj^^ a ^ a and hj^ )i ^ a ^ a from now on. 
After the equations of motion have been obtained from 
the Hamiltonian, one may identify primed and unprimed 
variables (e.g., the objects 1 and 1'), which in general re- 
quires another application of regularization techniques. 

The formulas for the radiation-reaction Hamiltonians 
Eqs. (3.34) and (3.35) can be simplified further. First, 
Eq. (3.34) may be written as: 



rrint _ 1 1.TT t 



r/ 3 T R STF 



(3.37) 



where the x-independent hj^jij was pulled in front of 
the integral and -B(4)ij is contracted with the symmet- 
ric trace- free hj^-. As explained previously, tiESy must 
be replaced by hj^ ji - a ^ a . The remaining integral in Eq. 
(3.37) is identical up to a prefactor to the definition of 
XU)iji cr - (3.28). Finally we obtain, inserting Eq. (3.18), 



rrint 
-"2.5PN 



16G 



X(4)ij X(4)ij 



(3.38) 



which is a well-known result (see [2] and references 
therein). The problem was reduced to the calculation 
of X(4)ij via (3.28). Remember that X(4)ij m this Hamil- 
tonian is explicitly time-dependent. 

We proceed with a simplification of the individual 
parts of Eq. (3.35). Analogous to the simplification of 
^25PN gi yen m the last paragraph we have 



1 



16ttG 
1 

16ttG 



16G 



xSi'Hui, (3.39) 



d 3 xh 



i,TT a^a'-iri ± 
(5)ij V (3)9(2), j 



■a-va tt 

16G ' 4 ' ij 



X(4)ij 1L 3ij 

(3.40) 



where Eqs. (3.29) and (3.31) were used. We may further 
write 

1 ' ,)3, lTT o^a'tTT a-j-q'qj matter _ g -a-ya' tt 



with the definition 



tt _ 4G / ,3 jTT a-s-a'-T/matter 



(3.41) 
(3.42) 
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The notation was chosen because of the similarity 
to n^i-,, cf. Equation (3.30). If the self-interaction con- 
tributions to the integral in (3.30) vanish, then E^y can 
be obtained from Ii-2ij by a relabeling of objects only. 
For the spin-dependent part of n 2 jj , this will turn out to 
be possible. The integral over B^jhJ^ in Eq. (3.35) 
splits into the following five parts, cf. Equation (3.19), 

1 



16ttG 
1 

16ttG 
1 

16ttG 
1 

16?rG 
1 

16vrG 



d 3 xfl a ^ 




TT a ~ > a 

~ 16G Uj 


X(4)ij , 


(3.43) 




a' n 


TT a— ^ a 

~ 16G 2<j 


X(4)ij , 


(3.44) 




a' i) 

U(4)ij 


TT a_ ~ *" a 


X(4)ij 7 


(3.45) 


d 3 xii a 4 ^ 


a' R 


= (<T) 3 (#- 


f 


(3.46) 


d 3 xQ a ^ 


a' TD 

B (i)ij 




+ 


(3.47) 



This agrees with [2] (with misprints corrected in [11]). 
It should be noted that no time derivatives are present 
in Eq. (3.3), so all time derivatives in Eqs. (3.38) and 
(3.54) are introduced by above insertions. Indeed, all 
these time derivatives should be understood as abbrevi- 
ations and be performed before the equations of motions 
are derived from the Hamiltonians. However, for time 
derivatives of primed variables it is irrelevant at which 
stage they are eliminated (These are actually all time 
derivatives except the one acting on Oij). One should be 
aware that an insertion of equations of motion leads to a 
recombination of PN orders, e.g., inserting the 1PN con- 
servative part of the equations of motion leads to 3.5PN 
contributions from H^ PN , cf. Equation (3.38). Further, 
one should notice that n 2 ij, R' , R" , Q'^, and Q'-j depend 
on both primed and unprimed variables by virtue of their 
definitions. 



Notice that here Tlaj, I^y, an d ^-3ij are independent of 
x. The relations Eqs. (3.33) and (3.18) were used in the 
last integral. The last two integrals were each split into 
two parts using Eq. (3.12) and the following definitions: 



R' 

R" = - 

Q'ij 



d 3 xT(±)ijH2~j a , 



1 



d 3 X 0( 2 ),i<£(2) jUSy" 



128ttG . 

'4, 

Q" _ l.l'... . A J I /./ » -"' 



A^nT rTT/ci ia-Xj' 
d X T(4)kl<>ij 0(2) 



256ttG 



to 1 1 r i ± kj 1 a— y, 

d X<p{2),k<P{2),l0ij 9(2) 



(3.48) 
(3.49) 
(3.50) 
(3.51) 



The fact that the explicit time derivative df x only acts 
on primed variables was used in (3.46) to pull it in front 
of the whole expression. Finally, it holds that 



1 



d 3 xh7Z^ a '^ 



16-G / 
with the definition 

° 13 = 16^G 



(5)matter ~~ X(A)ij J 

(3.52) 



(5) matter 



(3.53) 



Summing up the contributions from Eqs. (3.39) - 
(3.41), (3.43) - (3.47), and the total time derivative of 
Eq. (3.52), one gets: 



rrint u 

^3.5PN " 



+ x^(ni« + n 2tf )-x^5'n3« 
+ xt^'(Qi ] +Q" 3 ) + (dn 3 (R , + R'') 



d r 
dt 



X(4)ij Oij 



(3.54) 



IV. CALCULATION OF THE HAMILTONIANS 



Up to formal 3.5PN order, the interaction Hamiltonian 
is given by Eqs. (3.38) and (3.54). The quantities enter- 
ing these expressions must be calculated by solving the 
integrals appearing in their definitions [see Eqs. (3.17), 
(3.28-3.32), (3.42), (3.48-3.51], and (3.53)). The leading- 
order source terms in the pole-dipole case entering these 
integrals read 



a 



PaiPajSa + P a (iS a (f))(k)dk5a 



#(4)ij - -87TG7(4)ij - g</>(2),i0(2),j 



(4.1) 

(4.2) 
(4.3) 



see [1] for more details. Here, m a (a — 1,2...) are the 
masses and 5 a — <5(x — z n ). (j>(o\ is proportional to the 
Newtonian potential of point-masses, namely: 



(2) = — 167TGA = 4G ^ 



r„. 



(4.4) 



where r a = |x — z a \. Notice that 0(2) is independent of 
the spins. The expression for 5(6)^ was derived in [1]: 
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P 2 5 P 2 1 

^PaiPajSa + ^^PaiP a j(p(2)Sa + ^^P a iS a (j)(k)S a ,k ~ J^P a lPajPakS a (l)(i)S a ,k 



4m£ 



4m£ 



5 

+ fy^ PaiSa U)(k) (<P(2)5a) tk + -^P a iS a (k){j)<l>{1),k$a ~ ^-PakS a ( k )(i)<t>(2),j$a 



1 



1 



+ 2 S <k)(i) ( V (3),k + V ?3),j) 6 » 
+ ~<Pl(i)<P(2),ij + l<t>2(4)4>(2)^ + ^<f>(2)<f>(2),i<f>(2),j + 2ffJJ ~ %fe) + 2jt (3),k V m + ^(3)* 



(4.5) 



(3),fc 



The field quantities entering Eq. (4.5) are equal to: 



been computed in [14] and reads: 



0i(4) = 



PaiS a (i)(j) ( 1 



m„r„ 



m a m b 



^(4)— 2G 2 ££ — 

Pa, 



'(3) 



= «E 



5, 



o»C7) 



P 1 

9 — P 4- << 



^(3) 



«E 



2P Qi 



2P„ 



, 1 \ 1 

VijPak [ J ~^PakJ 'a,ijk 



(4.6a) 
(4.6b) 

(4.6c) 



,3 - 

(4.6d) 



(4.6e) 



'a(fc)(i 



a / ,fc 



1 



— -5, 



a/ ,kj 



o.(k)(j) 



a / ,ki 



where r ao = |z a — z \. Notice that for nonspinning sys- 
tems the result in [2] is reproduced. Further notice that 
nfy i does not depend on spin. Finally, the spin correc- 
tion to the field momentum is given by: 



T *3 

(5)mattcr 



4irG 



- ~ E ~^r P ^ P a(iSaU)Kk)Sa , (4.7) 



to the required order. 



A. Spin-dependent part of ftS 



The explicit solutions for the point-mass, i.e., spin- 
independent, contributions to hj^ can be found in [2, 
8, 11, 57] (but notice that [2] contains some misprints). 
The spin part of ^Suj> arising from the spin-dependent 



iTTspin _ n ST^ PanS a (k)(l) 



(45 k (i8j) n di - 2SijS k ndi)- 



(Skndidjdi - 25 k { i d j )d n di)r a 



(4.8) 



where we use the superscript "spin" to denote the spin- 
dependent part of a quantity from now on. In order to 
obtain the spin contributions to the radiation-reaction 
Hamiltonian up to formal 3.5PN order, we also need to 
compute the spin part of hj^j and hj^j ■ fy^y 1 "" ^ s eas Y 
to compute. From Eqs. (3.28), (4.3), and (4.2), we have 



spin 



SG \ ~s 



PaiS a (j)(k) 



d 3 xd k 5 a 



STF 



= 0, 



(4.9) 

and thus also h, 5 s^P m = [see Eq. (3.18)]. There is no 

spin contribution to the 2.5PN hj^ , which is the rea- 
son why the leading-order source terms (4.2) are not 
sufficient to derive the leading-order radiation-reaction 
Hamiltonian. hj^: spin would be more difficult to derive, 
but it is not needed in our calculation of the leading-order 
radiation-reaction Hamiltonian with spins, so we will not 
discuss it in the present paper. 

Analogous to Eq. (3.19), we decompose the solution 
for hj^ff m into several parts, 



"(7)»j 

TT spin 



,-rspin 

where the following definitions are used: 



n 



spin 



TTSpin 



,-rspin 
LL 2ij 



TTSpin 
LL 3ij 



n 



spin 



l 

5^ 
_4G 

5~ 

1 

~5n 
2G 



d 3 xBf£ Fspin . 

(6)iJ 



/j3„, lTT spin t, matter 

3 ( t t% spin / 
X { V (3) 0(2) ,0 



(4.10) 

(4.11) 
(4.12) 
(4.13) 



5 T J i} I dV|x-xf7X(x',t), (4.14) 



and obviously Ql 



0, cf. Equations (3.29) - (3.33) 



source terms in Eq. (4.2) via Eqs. (4.3) and (3.17), has and (4.3). These integrals yield the results: 



■prspin 
L1 lij 



yjspin 
LL 3ij 



4G 2 fir 

"rEE I p~ 3 ( naf) ' P b)"ab( n afc 5 'a( 4 )(fc) + Kb S a(j) (k) ) ~ 3P b fe (n^5 a (l) + n*(,S y) (fc)) 

- 3n^(P hj 5 a(l)(fc) + P b iS aU)(k) ) + 4(3< b n^ - 8 lj )n k ab P bl S a{k){l) +—^ r P a k{n J ab S a{l){k) + < 6 S y)(fc)) 

J m « r afc L 

+ (4^ - 6< b n J Qb )n^P ai S' Q(fc)(/) + 4n£ 6 (-P oi £ oW(fc ) + P ai S aU){k) ) - (3< b n^ b - %)5 b(fc)(z) 

+ 3n a6«6 5 '6( 4 )(0 + <6-%)(J)) + 3 (^i - 5 <6"ab) n afc<6 S '6(n)(0 

2 -Pafe(raa 6 5 a (j)( fe ) + n J ab S a (i)( k )) + n ab (P a iS a (j)( k ) + PajS a (i)(k)) 



4G 2 m b 1 f 

K 2-^2.^ ma r 2 b 



yjspin 

u 4« - 15 



+ 3(n Qb • P )n* 6 (n* 6 5 y) (fc) + n ] ab S a{i){k) ) + (<% + Zn l ab n 3 ab )n k ab P al S a{k){l) 

4G 2 \ \ m b u a in \ 

-^rl^l^ — n ab{n ab S a(i)(k) + n ab S aU)(k) ) , 
m 



^ — [-Pafe(«i»S'a(i)(fc) + ^a^a(i)(/s)) ~ 2 n a(PajSa(i)(k) + PaiS a ( 3 )(k) + 5ijPalS a ( k )(l)) 

lit* n L 



(4.15) 
(4.16) 



(4.17) 
(4.18) 



where n a = (x - i a )/r a and n Qb = (z Q - z b )/r ab . Notice 
that it holds 



T-rspin 

ll Uj 



T-rspm 
' U 2ij 



nspm _ 
L1 3ij ~ 



4G 

rspin 

~~5 *i 



(4.19) 



at the considered PN order, where 1^ is a multipole mo- 



I 

spin-orbit and spin(l)-spin(2) equations of motion [28- 
30], provided in this paper by Eqs. (5.2) and (5.3) later 
on]. 

The spin part of the formal 3.5PN order interaction 
Hamiltonian Eq. (3.54) can be written as: 



ment of the far-zone expansion of and can be ex- fj 
pressed as a double time derivative of a very compact 
expression, see Eqs. (6.15) and (6.18) in [1]. 



B. Derivation of spin contributions to 2.5PN and 
3.5PN interaction Hamiltonians 



rint spin 

3.5PN — 



16G 



K ll uj + l hij ) X(4)ij 

, s.a-Ki' ("rrspin , fj-spinx -a-Xi'TTSpin 
+ X{A)ij \. LL lij + vv 2ij I ~ X(4)ij L hij 

+ X%iQT n + (d° x ) 3 (R' spiD + i?" spin ) 

(4.22) 



d_ 
dt 



X(A)ij U ij 



When taking into account the fact that Eq. (4.9) tells 
us that X(4)ij = 0j we immediately see that the formal 
2.5PN order interaction Hamiltonian Eq. (3.38), 

5 . „^„' 



tt mt 
-"2.5PN 



^QQ X{A)ij X(A)ij 

has only the well-known point-mass contribution [57 

" 2 



(4.20) 



4G^ 
X(4,)ij - -jig 2^ 



"(Pa^ij SPaiPaj) 



GE— (V3«> J J 



Tab 



(4.21) 



but no direct spin contribution. However, indirect spin- 
contributions arise from Eq. (4.20) via the time deriva- 
tive therein and first appear at the formal 3.5PN level 
[after taking into account the leading-order conservative 



where we used X(4)ij = an d Qij Pm — 0- The latter 
is trivial from Eq. (3.51), as only the spin-independent 
potential </>( 2 ) appears there. n^ n , nS", n^ 11 , and 
were already derived in Sec. IV A. The missing 
quantities fi^f, Q[f in , i?' s P in , i?" s P in , and 0*f m can be 



jspin ^yspin 

obtained from 



rrspin 
LL 2ij 



Q/spin 
ij 

j^fspin 
j^ffspin 



d 3 xhJLT na ^ a 'n? lllttei 



"(4)ij 



4G 
~~5~ 

1 f ^3^,7-spin rTTfei ,a-)-a' 
4 '{4)kl d ii ( P{2) 



2 , d\{T {i)l3 ni-; a 'r^ 



l (2) 



1 1 

8 16ttG 



spina— > a 



j3 ij spin 
16ttG / (5)matter ' 



(4.23) 

(4.24) 
(4.25) 
(4.26) 
(4.27) 
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using Eqs. (3.42), (3.48) - (3.51), and (3.53). We also 
split R' s P in into three parts, 

R ,spin = rfspin + R ,n P in + R ,spm ^ (4 2g) 

i?'i Spin = -\j d " x ^ n ™ a "°' . ( 4 - 2 9) 

R'r n = -\J ^™H~" ' , (4.30) 



p/spin _ x I ,3 ^piiiT-rspina-s-a' 



(4.31) 



Here, PM denotes the point-mass parts of a function. 
The quantities entering above integrals will be all given 
in the present paper, except for IT™, which can be read 
from Eq. (36) in [2] using Rij = dfll^ij. The results of 
the above integrations read as follows: 



T-rSpin 
LL 2ij 



Q 



/spin 



/spin 



p/spin 



p/spin 
^3 



R 



AG 2 x m a 1 



E 2p a'k{n l aa ,S a , (j){k) + n 3 aa ,S a , {i){k) ) - n k aa ,{P a ,iS a , {j){k) + P a 'j 



3(n aa / • Pa')«L'«a' S 'a'0)(fc) + n L' S a' (i)(k)) ~ ( S U + 3 Ka'Ka') n aa' P a'lS a '(k)(l) 



G x - r?y 1 

4 ^— m„ 



2-Pafe« a '<5'a(j)(fe) + n i a > S a (i) (k)) ~ n aa> ( P aiS a ( 3 ) (k) + PajS a (i)(k)) 



3(rw • P a )n k aa ,{n l aa ,S a{j)(k) + n J aa ,S a(i)(k) ) - (*y + 3< Q ,n^ a ,)«aa- p aiS'a(fe)(;) 



G 
15 



E^ 



(00') 



m a >m a 



G \ - m a ,m b > 



E 



7 " m 

6r?, 



-P 'a' n a> a P aj (j^-a' a ' P a')Pa' iPaj 2(P a ' • P a )n a , Q P a 'j 

2r„> 



17n a » b /P a j 



07* , r 

-/a <, b ,P aj +2(n a , a .P a )<, fc ,<, a 



17(iv 6 , • P a )n l a , b ,n 3 a , a + 7n a , a l P a] 
8r a > a 



S a ' (i 



(00) 



//spin 



' a'b 

4G r aa / 
15 ^— ' m a >m a 

a, a' 

— ^ 77 l a 777, 5 T a ' a 



a'6' 



t- 2 P — r z n l P 

1 a'a^a'a 1 aj 'b'a"'a'a 1 aj 



2 „i 



~P a n aa' Pa' j 2(P a ' ■ P a ) n a a' Paj + ( n aa' ' Pa)-fa'i-f(ij 



WO') 



-^Pa' kPaiSa' (k)(j) ~ (Pa' ' Pa)<Sa' (i) (j) ~ Pa' iPakS a ' (k) (j) 



EE 



15 ^.t? m a' r ab 



Sa' 



(i)U) 



T-a'a'Pa'j ~ 2(n ab • P a ')K'a n ab ~ ( U a'a ' na^riafc'Pa 



(4.32) 



(4.33) 



(4.34) 



(4.35) 
(4.36) 

(4.37) 



sp in = £ 



-^2Pak{PaiS a {k)(j) + Pa 3 S a {k)(i)) ■ (4.38) 



The term in Eq. (4.34) containing 17n l a , b ,P a j actually 
cancels if the sums over a' and b' are performed and may 
therefore be dropped. 

Notice that n^", n^", and Q*? m are given by almost 
identical expressions, cf. Equations (4.16), (4.32), and 
(4.33). This is not accidental, but due to similarities of 
their defining integrals. With the source mass density 
given by Eq. (4.1), we obtain from Eq. (4.12): 



T-rspm 



4Gv^ 



, TT spin 



Similarly, Eq. (4.23) leads to 



T-rspin 
LL 2ij 



4G 



E/.TT spins. 
m (/i (4)i / ) 



(4.39) 



(4.40) 



Notice that in this expression no regularization is needed 
for taking x = z a , as primed and unprimed objects are 
not identified yet. In contrast to that there may be con- 
tributions from Hadamard regularization in Eq. (4.39). 
However, for the spin-dependent part, no such contri- 
butions appear (in contrast to the nonspinning case in 
[2]), which explains the great similarity between n^y" 

and n^j n . Further, insertion of Eq. (4.4) into Eq. (4.24) 
leads to: 



Q 



/spin 



d 3 xhJl. spin (n^ ttcI ) a - 



1 \ " j, TT spin 



"(4)ij 

h 



(4)i? 



(4.41) 



after performing several partial integrations and using 
Eqs. (4.3) and (3.17). Here also no regularization is 
needed. The similarity to Eqs. (4.23) or (4.40) is ob- 
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vious. The difference is simply an overall factor and a 
mutual exchange of primed and unprimed variables. 

V. ENERGY LOSS OF A BINARY SYSTEM 

A. Derivation of the energy loss from the 
Hamiltonian 

The instantaneous (near-zone) energy loss of a two- 
body system due to gravitational radiation can be writ- 
ten in the form (see, e.g., [2, 11]): 

/•inst Qox/rrint . rrint \ fK -l\ 

^S.SPN — °t \ n 2.5PN + -"3.5PNJ ■ K - 1 ) 

Notice that this energy loss is gauge-dependent in con- 
trast to the energy flux at infinity. 



We substitute Eqs. (4.20) and (4.22) into Eq. (5.1) (For 
the point-mass part of H™^ PN this was already done in 
[2]). After that, we need to eliminate the time derivatives 
in Eq. (5.1) using the leading-order spin-orbit, spin(l)- 
spin(2), and Newtonian equations of motion derived from 
the corresponding Hamiltonians [see, e.g., Eqs. (7.28) and 
(7.29) in [14]), 

z\, = ^ - ^;^F (3m 2 ,S v( j )(i) +im v S 2 >(j)({)) , 

(5.2a) 

4 = (1^2), (5.2b) 



Pi' 



Gmym 2 ' j Gn\i 2 i j 9 m 2 ' 



V2' 



15 

n>2> 

9 my 
2 ra 2 ' 



(ny 2 > • Sy) (n V2 > ■ S 2 
(p 2 ' x S 2 ') • rii-2 



(Pi' x Sy) • n V2 
3 



ry 2 < 



(Si/-S 2 /)+6 



(P2' x Si-) • nx'2' 
(pi' x S 2 ') • ni' 2 . 



G J 3m 2 ' 
2 my 

3 mi 



( P y X S!') 4 +2(p 2 ' X Si 



ri>2> 



(ni' 2 '-S 2 ')Si'-2(pi' xS 2 ' 



( P 2'XS2')' (ni' 2 '-Si')^, , 

2 m 2 ' ri'2' 



(5.3a) 



p\, =(1^2). (5.3b) 

Note that because the 2.5PN order Hamiltonian does 
not have spin contributions as we discussed in Sec. IV, 
we do not include the 1PN point-mass terms because 
substituting them into the 2.5PN Hamiltonian only pro- 
duces point-mass terms at 3.5PN order, while substi- 
tuting them into the 3.5PN Hamiltonian only produces 
4.5PN terms which is beyond the scope of this paper. 

At this point, we no longer need to distinguish the dif- 
ference between the primed and unprimed variables. Us- 
ing the Hadamard regularization method, we remove the 
singularities produced by the limit iy — > Zj and z 2 > — > z 2 
and obtain an expression of the energy loss in terms of 
Zi( 2 ) and Pi( 2 )- By realizing that z a = v a , we may use 



G n J 

p\ = m\v\ - (3m 2 Siy)(i) + 4mi Sat/Ho) > ( 5 - 4a ) 



P\ = (1 ^ 2) , (5.4b) 

to express the particle momenta p a in terms of the parti- 
cle coordinate velocities v Q , which can be easily obtained 
from Eq. (5.2). Here, r — ri 2 , and n = rii 2 . Note we do 
not include the 1PN point-mass terms in this expression 
for the reason described above. 

To put the energy loss into a more convenient form, 
we rewrite the individual masses mi , m 2 into the total 
mass of the system M = mi + m 2 , the reduced mass 
/i = mym 2 jM^ and the symmetric mass-ratio parameter 
j] = fj,/M using the relations (assuming mi > m 2 ): 

mi = £ (l + y/1 - 4ry) , (5.5a) 

m 2 = ii (l - y/1 - 4ry) . (5.5b) 

We also transform the individual coordinate velocities 
of each particle into the center of mass frame using the 
relations: 
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Vl = 



v 2 



^ r 

1 + v/l - 4t] 4r 2 

-2??v G r 

1 - y/1 - 4rj At 2 



(nx Si)(-1 
(nx Si) ( 1 



\/l - 4t] J + (n x S a ) (l + y/l - ^ 
V 1 - 477) + (n x S 2 ) (l + V 1 - 4 ?/) 



(5.6a) 
(5.6b) 



where v = Vi — v 2 is the relative velocity, S a is the indi- 
vidual spin. Notice that here we do not include the 1PN 
point-mass terms (see, e.g., Eq. (3.13) in [11]) because 
the 1PN corrections of v a can only produce 3.5PN terms 
in the flux when substituted into jC<2*5pn' which is in- 
dependent of spins, therefore the 1PN point-mass terms 



r 



in v a do not contribute any spin-dependent terms at the 
formal 3.5PN order. 

After eliminating the coordinate velocity vi and v 2 
by means of Eq. (5.6), the spin-orbit and spin(l)-spin(2) 
(S1S2) part of the instantaneous energy loss £<3* 5PN can 
be written as: 



„ inst 4 G W G 2 M 2 o . GM , . , 2 2 J 

AI'spn = +45(n • v) 4 - 60(n • v)V + llv* + (-9(n • v) 2 + llv 2 ) , 



-.inst.SO 
-<3.5PN 



15 r 3 
G 2 M 2 n 2 



(5.7a) 



15r 5 



(Ln • i) 



(Ln • S) 



140- 



G 2 M 2 

74 — + 420(n • v) 4 - 510(n ■ v) V + 66v 

G 2 M 2 



GM 



GM 



(54(n • v) 2 + 22i 



+ 840(n • v) 4 - 840(n • v) V + 96v 4 + (336(n • v) 2 + 84v 2 ) 



(5.7b) 



C 



inst , SlS2 _ 2 G 2 Mr\ 



<3.5PN 



15 r 5 



(Si-S 2 ) 



_ 12Q(n . v) 2 w 2 + 24u 4 + GM_ / i92v2 _ ^ ^ 



(Si • v)(S 2 • v) 



184— - 450(n • v) 2 + 138v 2 
r 



[(n • v)(n • S 2 )(Si • v) + (n • v)(n • Si)(S 2 • v) 
G 2 M 2 



C 1 A J 

-546 — - + 1785(n • v) 2 - 1005w 2 



(5.7c) 



(n- Si)(n- S 2 ) 



- 36- 



5670(n • v) 4 + 4620(n • v) V - 39(V 



+ ^(-600« 2 + 1536(n-v) 2 ) 



with v = |v|, S = Si +S 2 , £ = (w2/mi)Si + (mi/rn 2 )S 2 
are the spin variables, and L^v = rn x v is the Newtonian 
orbital angular momentum per reduced mass. 



B. Comparison with other results 

References [25, 26] recently computed, using the 
method of direct integration of the relaxed Einstein equa- 
tions [59, 60], the leading-order spin-orbit and spin(l)- 
spin(2) equations of motion and the corresponding energy 
loss in harmonic coordinates. In this subsection, we shall 
prove that our result is actually equivalent to the results 



in [25, 26]. 

In order to compare the instantaneous energy loss, we 
first need to find the transformation between our ADM 
canonical variables (z a ,v a = z a , S a ) and the "harmonic 
coordinate" variables (y a ,V a = y a ,S^ w ). Because the 
quantity we are comparing is the energy loss £<35pn 
at formal 3.5PN order, which is only one formal order 
higher than the leading-order energy loss £<1*5pn caused 
by the quadrupole radiation of point-masses, the coordi- 
nate transformation we are looking for only needs to be 
accurate up to formal 1PN order. 

It is well known that for the point-mass case the ADM 
coordinates are equivalent to the harmonic coordinates 
at 1PN order in that they result in identical equations of 
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motion. In addition, the spin-dependent part of the for- 
mal 1PN accurate transformation z a (y Q , V a , S^ w ) can 
be derived from the well-known transformation between 
different spin supplementary conditions (SSC) (for de- 
tails, see, e.g., [25]). Namely, for a specific SSC param- 
eter k, which is used to fix the center of mass of the 
particle, we impose the condition: 



s: 



0. 



(5.8) 



where k typically has the value 1, 1/2, or 0. The relation 
between the center of mass for each value of k can be 
written as: 

(4) (fe0 = (4) (fc) + ^ #Vj (fe) - (5.9) 



It is straightforward to show that at formal 1PN order 
the SSC in our calculation leads to k = 1/2, which is 
identical to the one used in references [25, 26]. Therefore, 
we have: 



z a (y a ,V a ,SWW )=yaj 
v a (y a ,V a ,Sr W )^i a = V a 



(5.10) 
(5.11) 



Reference [37] has shown that the difference between 
the spin parameters S a used in the ADM formalism and 
the ones used in the harmonic coordinates calculations is 



of formal 2PN order. In other words, the transformation 



Sa(y Q ,V a ,sww ) = sr 



j WW 



(5.12) 



can be used in this paper. 

From Eqs. (5.10) - (5.12) we know that our ADM 
canonical variables are actually equivalent to the har- 
monic gauge ones at the considered PN order. Now we 
are not comparing with the harmonic gauge energy loss 
given in [25, 26, 60] directly, but with the far-zone energy 
flux, which was shown to agree with the former (up to an 
nonphysical total time derivative). When comparing our 
result Eqs. (5.7) to the far-zone flux [£<3.5pn] ar " zonc , for 
the purpose of this paper, only the parts 



<3.5PNj 



[C 



<3.5PNj 



i far-zone 



[£<2.5PN + £<3.5PN 
i /»SiS2 "linst 

' ^-o.spnJ ) 

[£<2.5PN + £<3.5PN 



C 



SiS 2 
<3.5PNJ 



1 far-zone 



(5.13a) 
(5.13b) 



are relevant to this paper, where for the instantaneous 
energy loss in ADM coordinates we substitute Eq. (5.7) 
and for the far-zone flux we substitute the expressions 
computed in [19] in harmonic gauge, 



<2.5PNj 



lf ar-zo„e 8 G A 



4^2 



15 



[(-ll(n-v) 2 + 12v 2 ] , 
(5.14a) 



|/-<3.5PN 



far- 



8 G 3 MV 
15 




Lm • S 



GM 



+ 18(n • v) 2 - 19v 2 



-12— + 27(n • v) 2 - 37« 2 
r 



(5.14b) 



C 



<3.5PN 



far-zone 4 (J° M T) 



15 r 6 



-171(n • v)(n • S 2 )(Si • v) - 171(n • v)(n • S X )(S 2 • v) + 71(Si • v)(S 2 • v) 



(n-Si)(n-S 2 ) 807(n-v) 2 - 504u 2 + (Si • S 2 ) [-165(n ■ v) 2 + Ulv 2 } 



(5.14c) 



It should be noted that the sources (on the right-hand 
side of these equations) are evaluated at the retarded 
time with respect to the flux (on the left-hand side), 
which is not explicitly denoted here. In contrast to the 
instantaneous near-zone energy loss, these results are ac- 
tually gauge-independent at the considered PN order, 
i.e., one gets exactly the same result from Eq. (6.22) 
in [1] within the ADM gauge. We already showed this 

in [1] for [£<3 5PN ] ar ZOn °, and we confirmed this for 

c r ~ 1 far-zone 

1 tar-zone n 

and 



<2.5PNJ 



S1S2 
<3.5PN 



too. 



It has been shown in [11] that the spin- independent 
result [£<2.5PN] mbt + [^3.5pn] an d the spin-orbit part 
of formal 3.5PN order [£f ( gp N ] 1118 agree with the results 
computed in harmonic coordinates up to a total time 
derivative, which is a pure gauge effect and vanishes after 
orbital average (see, e.g, [27] and [61]). 

For spin-dependent instantaneous energy loss, it is pos- 
sible to write the difference between Eqs. (5.13a) and 
(5.13b) as a total time derivative using the identities in 
Appendix A, which has already been presented in Ap- 
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pendix F of [25] and Appendix A of [26]. Taking into 
account the leading-order point-mass and spin contribu- 
tions it holds that 



£rSO i /»SiS 

<2.5PN + *-<3.5PN + L - 



<3.5PN 



far-zone 



d_ 
dt 



-B2.5PN + -B3.5PN + ^3.5PN 



, (5-15) 



where 



£<2.5PN + £<3.5PN 



/*Si S2 
^S.SPN 



-C/2.5PN = 5 (n • v) 



44 , 12, N2 8 GM 

— v A (n-v) , 

15 5 v ' 15 r J 

(5.16a) 



El 



SO 

3.5PN 



E: 



3.5PN 



G 2 M 2 if 



(n-v) 



(L N • S) 



28 GM 



( n . v )2_f2 2 



15 



, 2 GM t . , 2 22 2 
(Ln-0( +4(n-v) 2 - — « 2 



G?7 



(n-S 1 )(n-S 2 )(n-v 



GM 



84 (n .v)^-52^- 44GM> ) 
or/ 



GM ( „ . . „ wo A / . , 2 38 2 22 GM 

— f (n ■ Si)(S 2 ■ v) + (n • S 2 )(Si ■ v) j f -22(n ■ v) 2 + — v 2 + — — 

16GM /0 .,_ w , 16 GM 2, v /0 
+ (Si • v)(S 2 ■ v)(n • v) + — v 2 (n • v)(S x • S 2 ) 



5 r 



5 r 



(5.16b) 



(5.16c) 



Note that even though the energy loss at 2.5PN order 
is spin- independent, it does need to be taken into ac- 
count when comparing the spin- dependent energy losses 
because of the spin-dependent terms in Eqs. (Ala) and 
(Alb), which are of formal 1PN order. 

It should be noted that the total time derivative on 
the right-hand side of Eq. (5.15) vanishes to the order 
in question when averaged over time. This means that 
the time average of near-zone energy loss and far-zone 
energy flux agree. Equations (5.16a) - (5.16c) should 
thus be interpreted as (gauge-dependent) energies that 
temporarily leave the near-zone, but never reach the far- 
zone and instead move back into the near-zone at a later 
time. Therefore, they have in average no effect on the 
near-zone energy loss. 

VI. CONCLUSIONS AND OUTLOOK 

Based on developments in [1] the leading-order PN 
spin-orbit and spin(l)-spin(2) radiation-reaction Hamil- 
tonians were calculated. Corresponding equations of mo- 
tion were already derived for the binary case in [25-27]. 
The Hamiltonians given in the present paper are even 
valid for arbitrary many spinning compact objects and 
present the dynamics in a compact form. The derivation 
was performed within the ADM canonical formalism [53] , 
which was extended from point-masses to linear order in 
the single spin of the objects in [1, 12-14]. The calcu- 
lation of the needed integrals and their regularization is 
analogous to calculations for nonspinning objects within 



the ADM formalism (see, e.g., [2, 8, 62]). In particular, 
we applied the Hadamard finite part and Riesz-formula 
based regularizations in the present paper (for the latter 
see also [63]). Some integrals were checked using Riesz 
kernels in arbitrary dimensions (see also [64]). 

The leading-order spin-orbit and spin(l)-spin(2) en- 
ergy loss was computed in the present paper from the 
explicit time derivative of the interaction Hamiltonian. 
This was compared to well-known results for the corre- 
sponding energy flux [19] as a check [In [1], the leading- 
order spin-orbit energy flux was rederived from the wave 
equation (3.9)]. This also proofs agreement with the en- 
ergy loss obtained in the harmonic gauge [25, 26, 60] and 
thus provides an important check of the ADM canon- 
ical formalism for spinning objects, which was derived 
only very recently [1, 12-14]. Notice that the interac- 
tion Hamiltonian in the form of Eq. (3.3) also gives es- 
sential contributions to the next-to-next-to-leading-order 
conservative Hamiltonians [48, 49]. 

The spin-orbit radiation-reaction Hamiltonian derived 
in the present paper, which is at 3.5PN when counted in a 
formal way, is actually of the order 4PN for maximally ro- 
tating objects (see also Appendix A of [1]). A derivation 
of all spin-dependent 4PN Hamiltonians for maximally 
rotating objects should be envisaged in the future. The 
most complicated Hamiltonian at this level is the conser- 
vative next-to-next-to- leading-order spin(l)-spin(2) one, 
but it has already been derived very recently (see [49], 
and also [50] for a corresponding potential). Notice that 
all Hamiltonians for maximally rotating black holes are 
known to 3.5PN order [48]. 
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Further, the leading-order spin-orbit and spin(l)- 
spin(2) radiation-reaction equations of motion can be ob- 
tained from the Hamiltonians derived in the present pa- 
per and compared with the results from [25-27] in the fu- 
ture. Primed and unprimed variables must be identified 
in the equations of motion, which requires further ap- 
plication of regularization techniques. Finally, one may 
transform the general equations of motion into secular 
equations of motion for the orbital elements, which has 
already been derived in [22-24] using energy and angular 
momentum balance. 
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The identities for total time derivatives needed to com- 
pare the instantaneous near-zone energy loss and the far- 
zone energy flux were already provided in Appendix F of 
[25] and Appendix A of [26]. They read 1 : 



d fv 2s r p 
dt V ri 



d ( v 2s r p . 



dt V ri 



L N = 



d ( v 2s f p 



dt \ ri 

d (v 2s r p 



dt \ r<i 
d ( v 2s r p 



.1^3 



dt V r<i 



r q- 



— < pv 4 — (p + q)v 2 f 2 — 2sr 



GM 2 GM p Gv 2 f 



pv 

r r 



2 r 3 



L N • (45 + 3£) 



Gr 



6s^ [r (Si • S 2 ) + (v • Si) (n • S 2 ) + (v ■ Si) (n • S 2 ) - 5r (n- Si) (n ■ S 2 ) 

3^[(Si-S 2 )-3(n-S 1 )(n-S 2 )] 1, 
fir 6 



r q+l 

pGv 2 ~ 



4 / \ 2-2 „ -2 GM 1 GM 

pv — (p + q)v r — 2sr pv 

r ' r 



f [ I^T-Ln- (45 + 30 J Ln - x ((v - ^rn) x (45 + 3f) 



6s 



Gr 
/ir 3 



f (Si • S 2 ) + (v ■ Si)(n • S 2 ) + (v • S 2 )(n • Si) - 5f(n • Si)(n • S 2 ) 



(Ala) 



(Alb) 



3p^ [(Si • S a ) - 3(n ■ Si)(n ■ S 2 )] L N - t( n x S i)( n ' S 2) + ( n x s 2)(n ■ Si) 

fir 



fir 

r q+l 
v 2s-2-p-l 

r q+l 
v 2s—2f,p—l 

r 9+l 



pv — (p + q)v r — 2sr 



pv — (p + q)v r — 2sr 



pv — (p + q)v r — 2sr 



. GM 2 GM 
pv 



r r 

. GM 9 GM 

pv z 

r r 

. GM 2 GM 

pv z 

r r 



x l x j + 2v 2 rrx^v 3 ^ 



2 • 

v l v j - 2GM V ^-x (i v j) \ , 

/ 1 • / 7 4 GM ■ ■ 

x v J + v tt v v J n n J 



(Ale) 
(Aid) 

(Ale) 



where r = (n • v), and s ,p ,q are non-negative integers. 



1 There was a misprint in Eq. (Al) of [26], we made appropriate 
changes in the expression here and marked its position by .... 
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